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C^ ■ Abstract 



O ■ We investigate lepton flavor violation in the context of intersecting D-brane models. We 

point out that these models have a source to generate flavor violation in the trilinear scalar 
^ couplings while the geometry of the construction leads to degenerate soft scalar masses for 

p ■ different generations (as in the minimal supergravity model) at the string scale. The trilinear 

ji. . scalar couplings are not proportional to the Yukawa couplings when the F-term of the U- 

D ! moduli contribution is non-zero. Consequently, the lepton flavor violating decay processes are 

'^ ' generated. Only other sources of flavor violations in this model are the Dirac neutrino Yukawa 

coupling and the Majorana couplings. The observed fermion mixings are realized from the 
^ ■ "almost rank 1" Yukawa matrices, which generate a simple texture for the trilinear scalar 

j^ . terms. We calculate the branching ratios of r ^ fi'j, yU — > 67 and the electric dipole moment 

of the electron in this model. We find that the observation of all the lepton flavor violating 
decay processes and the electric dipole moment will be able to sort out different flavor violating 
sources. 



1 Introduction 

The standard model is well established to describe physics below the weak scale while it also 
has a number of parameters, especially in the flavor sector. Indeed, the patterns of masses and 
mixings for quarks and leptons are not very simple and should be explained in a fundamental 
way. Thus, one expects that there exists more fundamental physics beyond the standard model 
and the masses and the mixings are described by some fundamental parameters. 

Supersymmetry (SUSY) is the most promising candidate of new physics. SUSY models can 
explain gauge hierarchy problems and suggest gauge unification such as SU{5) grand unified 
theory (GUT) with successful gauge coupling unification in the minimal extension of the SUSY 
standard model (MSSM). However, SUSY does not solve flavor puzzles. Rather, it increases the 
number of parameters with flavor indices to more than hundred in general. Nevertheless, people 
are not discouraged to consider SUSY models since the SUSY breaking parameters with flavor 
indices are constrained to suppress flavor changing neutral currents (FCNC) pp. Actually, one 
expects that the FCNC suppression may be realized by a flavor symmetry, which may give us 
a hint of the fundamental physics for flavors. 

The minimality of the SUSY breaking parameters is assumed in the minimal supergravity 
(mSUGRA) mediated SUSY breaking scenario [2 : SUSY breaking scalar masses are universal 
and the scalar trilinear couplings (A-terms) are proportional to the Yukawa couplings. The 
degeneracy of the SUSY breaking masses corresponds to the U{3)l x U{3)r flavor symmetry. 
On the other hand, it is hard to relate the proportionality of A-terms to such flavor symmetries 
since the Yukawa couplings themselves break the symmetries. 

The fundamental questions for the flavor sector are the following: 1) Why do fermions 
replicate with different masses? 2) Can we explain the pattern of the masses and the mix- 
ings for quarks and leptons? 3) Why does the flavor symmetry seem unbroken in the SUSY 
breaking mass terms, while the fermion masses break the flavor symmetry? 4) Is the A-term 
proportionality feasible? When is this proportionality feasible? How does it look like if it is 
not proportional? 

The intersecting D-brane models [3 El |3] may answer such questions. The string theory can 
describe the particle field theory as an effective theory, and thus, in principle, it has a potential 
to calculate all the parameters by using a few fundamental parameters. Indeed, the intersecting 
D-branes are interesting approaches to construct the standard model. The A^ stack of D-branes 



can form U{N) gauge fields as zero modes of open strings attaching on the D-branes. Open 
strings can be attached at the intersection between the N stack and the M stack of D-branes, 
and massless chiral fermions belonging to (A^, M) bi-fundamental representation can appear. 
Such a situation is very attractive to obtain quark and lepton fields not only in the standard 
model 6j but also in the unified models [ZlEl- When the extra dimensions are compactified by 
torus such as T^ = T'^ xT"^ xT'^ , the intersecting point of the D-branes can be multiplicated, and 
thus the fermions are replicated. The number of generation is therefore a topological number. 

In addition to the realization of the standard-like models, the effective supergravity La- 
grangian is calculable [HI CHI HH 1121 ■ The Yukawa coupling is obtained as an open string 
scattering for the triangle formed by the D-branes. The couplings are described as e~^^, where 
the triangle area A is formed by three intersecting points. For the toroidal compactification 
models, the Yukawa couplings are written as theta function of geometrical parameters including 
instanton effects. In simple models, the Yukawa matrices are written in the factorized form 
i/ij = xfx^ jni E] ■ This originates from a geometrical reason that the left- and the right-handed 
fermions are replicated at the intersecting points on different tori, and the Yukawa couplings 
are given as an exponential form of sum of the triangle areas. As a result of the factorized form 
of the Yukawa coupling, the Yukawa matrices are rank 1, and thus only the 3rd generation 
fermions are massive. In order to construct a realistic model, this issue of Yukawa matrices 
needs to be resolved and several possibilities have been considered [121 El . The Yukawa ma- 
trices can be hierarchical when the Yukawa matrices are "almost rank 1" by including higher 
order effects or quantum corrections. The exponential suppression of the Yukawa coupling is 
also available to obtain hierarchical masses. In Ref.^l], we have discussed that the observed 
patterns of fermion mixings can be easily reproduced if the Yukawa matrices are almost rank 1. 

The Kahler metrics of the zero modes can also be calculated as string scattering amplitudes 
[TU] in terms of the moduli fields: dilaton S, Kahler moduli T, and complex structure moduli 
U. The Kahler metrics are diagonal for the zero modes, and the metrics for the bifundamental 
fields are determined by the brane configuration parameters such as the relative angles of 
the D-branes. Since the relative angles are common when the fermions are replicated at the 
intersection of the D-branes, the Kahler metrics are flavor invisible. Consequently, the SUSY 
breaking scalar masses are same for different generations. So, the flavor symmetry of the scalar 
masses can originate from the brane geometry. 

When the Kahler metrics remain same for each generation, the Kahler connection parts 



(which are the derivatives of Kahler metric) of A-terms are common. Thus, the non-proportional 
part of the A-term is only the derivative of the Yukawa coupling. The Yukawa couplings which 
are given as theta functions depend only on the fZ-moduli ^2]; neither on the S nor the T- 
moduli. As a result, the trilinear scalar couplings are proportional to the Yukawa coupling when 
the F component of the [/-moduli is zero. If F^ ^ 0, the non-proportional part of A-term is 
acquired, which is proportional to the derivative of the Yukawa matrices. 

In this paper, we emphasize the degeneracy of the SUSY breaking mass terms and the 
[/-moduli contribution of the trilinear scalar couplings. These contributions are related to the 
flavor violation and we will study the lepton flavor violation (LFV) processes since the flavor 
violation in the lepton sector produce much more stringent constraint rather than in the quark 
sector. The LFV processes, such as /i — i> 67, r — > /i7 and r -^ 67, are not yet observed, but we 
have bounds on the branching ratios of these decay modes ^H] ■ The observation of these decay 
modes would provide information of flavor violation in new physics. It is pointed out that the 
processes are accessible for SUSY models [16^. In the mainstream of theoretical calculations, 
the flavor violation in the SUSY breaking parameters is assumed to be absent at the cutoff scale 
in the mSUGRA. The source of flavor violation originates from the Yukawa couplings such as 
the Dirac neutrino and the Majorana couplings. The SUSY breaking parameters at the weak 
scale can acquire flavor violation through renormalization group equations (RGE). Indeed, the 
off-diagonal elements of the left- and right-handed slepton mass matrices are generated, and 
the slepton-gaugino loop diagram provides the LFV processes. In the intersecting D-brane 
models, the flavor degeneracy of the slepton masses at the cutoff scale is realized naturally 
and the [/-moduli contribution of A-terms can be one of the major sources of flavor violation. 
We will calculate the branching ratios of different LFV decays and the electric dipole moment 
(EDM) of the electron, and study whether we can learn the origins of flavor violation from the 
forthcoming experiments. 

This paper is organized as follows: In section 2, we will study the low energy effective action 
of the zero modes for matter fields in intersecting D-brane models. In section 3, the realization 
of neutrino mixing angles are obtained in the context of "almost rank 1 Yukawa matrix" . In 
section 4, we investigate the sources of LFV. In section 5, we will calculate branching ratios of 
the LFV decays and the EDM of the electron, and compare the results with different setups of 
flavor violation. Section 6 is devoted to conclusions and discussions. 



2 Effective action in the intersecting D-brane models 

Our purpose is to deal with the flavor physics and study its phenomenological implications in 
the intersecting D-brane models. We will consider models with intersecting D6-branes in the 
type IIA theory [5 , which may be equivalent to the models with magnetized D9-branes and 
D5-branes in the type IIB theory. One can also apply our work to the D7-branes in the type 
IIB theory in which supersymmetry breaking soft terms can arise from 3-form fluxes J2]- In 
this section, we will explore the flavor sector of the models without concentrating on the details 
of any individual model. 

The MSSM-like models can be constructed easily by introducing three sets of intersecting 
D-branes. For example, in the type IIA orientifold models with T^ jZ^ x Z^ and with the 
intersecting D6-branes, the A^ stack of D-branes can form SU{N/2) gauge fields. Massless chiral 
fermions belong to the {N/2, M/2) bi-fundamental representation can appear at the intersection 
between the A^ stack and the M stack of D-branes. So, introducing a, b, c branes for f/(4)c, 
U(2)l and U{2)r respectively, we can obtain Pati-Salam-like model [TTj with quark and lepton 
fields IS]. The SU{A)c and SU{2)r symmetries are broken to SU{A)c ^ Sf/(3)c x U{1)b~l and 
SU{2)r —> U{1)r by brane splitting [TH], and U{1)b-l x U{1)r is broken down to f/(l)y by the 
Higgs mechanism. When the branes are parallel to the orientifolds, the USp gauge symmetry 
arises. The SU{2) gauge symmetry in the standard model can originate from the USp brane. 
In order to eliminate the RR tadpole, extra branes are often needed and they may form hidden 
sectors p!^ . 

Since the extra dimensions are compactified to T^ = T"^ x T"^ x T"^ , the D-branes intersects 
multiple times and the generations of the fermions are replicated. The intersection numbers 
lab for Da and Db branes are topological invariant and can be given by the wrapping numbers 
and the total magnetic fluxes (n^, m^), 
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lab = nK"^fe - KO. (1) 

r=l 

where r represents index of each torus. In a simple choice to obtain three generations, the 
left-handed matter \1/"'' and the right-handed matter \l/™ are often replicated on different tori. 
The important implication of the family replication is that the Kahler metric is flavor 
diagonal: 

K = k{M, M) + Kab{M, M)^f ^f + K,a{M, M)^™^f + . . . , (2) 




Figure 1: A sketch of brane intersections on a torus. 



where M stands for the 5, T and U moduh and the index i is for the flavor index. In addition 
to the flavor diagonal nature, the Kahler metric for ^""^ does not depend on flavor indices for 
lab = lea = 3. The Kahler metric Kab{M,M) is determined by the relative angles 6'^^ of the 
D-branes [131111112]: 



Kab oc e^^ n(^- + Ur 



\ 



r(l - Ur 



(3) 



where 04 is a 4-dimensional dilaton and Ur = Oab/'^, which are functions of S and T moduli. 
The moduli dependence of the Kahler metric is determined by the geometrical parameters, and 
thus the metric is flavor invisible. Therefore, the SUSY breaking scalar mass for the left-handed 
matter is given below ^20| and it has the flavor degeneracy: 



m 



ab 



m 



3/2 



Vo-Y.F^'F''dMdN\ogKab 



(4) 



KLN 



where m^i2 is a gravitino mass and Vq is the vacuum expectation value of the scalar potential. 
The scalar mass squared m1^ for the right-handed matter can be similarly written and can also 
have the flavor degeneracy. Since the relative angles for ab and ca can be different, the scalar 
masses are not necessarily universal for different representations of matter. We note that the 
flavor degeneracy can be broken when 2+1 decomposition of the generation is considered, for 
example, lab = 2 and lab' = 1, where b' is a orientifold reflection of the brane b. We will choose 
lab = Ica = ^ from uow ou and there is a ^(3)^^ x f/(3)ij flavor symmetry in the SUSY breaking 
scalar mass terms at the string scale. We emphasize that such a mSUGRA-type flavor structure 
can be obtained by a geometrical setup of the D-branes. 

The Yukawa coupling vl>'^^vl>^'=vl>™ is induced by the three-point open string scattering. When 
the left-handed matter \1'"* and the right-handed matter ^'^'^ are replicated in different tori, the 



Yukawa coupling is factorized: 

Y,,=x^{U{)xfm. (5) 

The x^ ' can be written by theta function [9^ with some geometrical parameters such as e 
shown in the Fig.l. Naively, these are given by e"*^"^ where A is the area of the triangle formed 
by the branes. The Yukawa couplings do not depend on the S and T moduli but depend on 
the complex structure moduli U. 

When the Yukawa coupling is factorized, the matrix is rank 1 and consequently U{2)l x 
U{2)f( flavor symmetry remains and the fermions of 1st and 2nd generation are massless. Surely 
such a situation is not viable, and there exists several discussion on this issue ^^. For example, 
we have suggested that the multi-point function of the string scattering including extra branes 
can increase the rank of the Yukawa matrix ^^. In this paper, we do not specify how to 
increase the rank, but we assume the factorizability of the Yukawa coupling at the leading 
order since this assumption leads to interesting phenomenological implications which we will 
see in the next section. 

The scalar trilinear coupling [A-term) is given as |2U| 

A, = F^ 



Km - Qm \og{KabKbcKca) ] Yij + dMYij , (6) 



and the coupling is proportional to the Yukawa coupling if F^ = 0. However, when F^ ^ 0, 
the flavor violation is generated, 

F'^^duJi, = F^'xfxf + F^^xfxf, (7) 

where x stands for the derivative by the U moduli. The U moduli contribution in the A-term 
can be the source of flavor violation. In this paper, we will emphasize the effect of the U moduli 
contribution. 

3 Application of "almost rank 1 Yukawa matrices" 

We will first discuss the mixing angles for neutrino oscillation in the context of "almost rank 1 
Yukawa matrix" [Tl] . 

The Yukawa matrix for the charged-leptons is given as Ye = Yo + 6Y where Yq is the rank 1 
matrix, and 6Y is a small correction to generate electron and muon masses. The rank 1 matrix 



can be expressed as 



Fn 



hi 



( C2 62 02 ) 



'^ C1C2 C162 Cia2 ^ 

61 C2 6162 ^102 
\ aiC2 0162 01^2 / 



(8) 



and Oj, fej, Cj can be given by theta function f^. Note that Oj, 6j, q can be rotated to be real by 
field redefinitions. 

Now let us work in a basis where the mass matrix of the light neutrino is diagonalized. 
Then the diagonalizing matrix of Y^ is the MNSP (Maki-Nakagawa-Sakata-Pontecorvo) matrix: 
^MNSP = Uf^, where UfYJJ'jl = Y^^'^^. The 3x3 unitary diagonalization matrix has three 
mixing angles, and those angles may be generically large since Oj, 6j, q are all order one pa- 
rameters. However, one of the three mixing angles of Ul(r) is unphysical in the limit 5Y —>■ 
since 1st and 2nd generation masses are equal to be zero and the U{2)l x U{2)ji flavor sym- 
metry remains unbroken. The small correction, 6Y, eliminates the degeneracy and the mixing 
angle of Ul{r) is fixed. Namely, the two mixing angles in Ul{r) are generically large and one 
mixing angle is determined by a small correction 6Y. For example, when the small correction 
is SY = diag(0,0,e), the Yukawa coupling becomes rank 2 and the eigenvector for the zero 
eigenvalues is oc (61, — ci, 0) for Ul- As a result, one can find that f/e3(= sin6'i3) is exactly zero 
in this example. The small correction needs to be more realistic to generate the electron mass, 
and then Ue3 can acquire a small non-zero value. Consequently, the two large mixings for solar 
and atmospheric neutrinos and the small mixing for ^13 is elegantly realized in this scheme. 

The approximate diagonalization matrix f/° is given as 



/ cos 6^ 



Ul 



sin6'f 



cos 9^ sin 9^ cos 9^ cos 9^ 





sin6'i 



\ 



(9) 



\ sin6'„sin6'^ sin6'„cos6'3 cos6'„ / 



where tan 6'^ = Ci/61 and tan 6*^ = Jhi + cf/ai. The right-handed U^^^ can be also described 
similarly. Then the MNSP matrix can be written as 



f/i 



MNSP 



vi*ul , 



(10) 



where V^ is a diagonalizing matrix of U^YgU^ = Yq '^^^ + Ul5YU'^ . 

In the quark sector, the CKM (Cabibbo-Kobayashi-Maskawa) matrix is written as Vckm = 
^L^L 5 where the unitary matrices are U^' Yu^JJ^ = Y^^^^. The Yukawa matrices are given 



as Y^' = xfxj + 6Yj^' . In a similar way in the charged lepton sector, the unitary matrices 
can be written in the form U^' = V^' U^. Since the left-hand part xf is common for both 
up- and down-type quarks, the large mixings in U^' get cancelled, and the CKM mixings are 
small: Vckm = ^l^l ■ Since the up-type quarks are more hierarchical than the down-type 
quarks, the CKM matrix is expected to be Vckm — '^l- If we have quark-lepton unification, 
we have a relation V^ ~ V^. Then the MNSP matrix is 

^MNSP ^ V^JkM^E- (11) 

This type of MNSP matrix is surveyed as an ansatz in Ref.|21j. 

Once the MNSP matrix is given in the form Eq. ()l()|l . we obtain the mixing angles for 
neutrino oscillation as follows |14j : 

sin 6*13 ~ sin 9^ sin 91^, 6'atm ^ 6*^, 6*301 - 0^ ± 9 13 cot 6'atm cos 6mnsp , (12) 

where 6*^2 is a mixing angle in VI , and 91^ and 92^ are neglected since they are expected to 
be small as in the quark sector. The atmospheric neutrino mixing is almost maximal and the 
solar mixing angle is large but not maximal since ci ^ bi <^ ai and tan 6'^ = Jbi + cl/ai, 
ta.n9g = ci/bi. 

The smallness of the neutrino masses are explained by the seesaw mechanism j221- We 
note that the large mixing angles between the charged-lepton and the neutrino Dirac Yukawa 
coupling can also get cancelled as in the quark sector. Thus, the favorable situation is when 
the SU{2)l triplet Majorana part is dominant in the type II seesaw 



flight = Ml- M^M^^M^^. (13) 

The Majorana couplings for both left- and right-handed leptons 

^hUAr. + i(/f e^e^A^- + /j^W^A?, + 72/^^=6= A],), (14) 

can be generated by multi-point functions in each torus |E]- 

4 Possible sources of lepton flavor violation 

In this section, we will describe the sources of the LFV processes such as ^ ^ e'-f and r -^ fi'y. 

In SUSY models, the LFV processes are described by loop diagrams. The charginos and 

neutralinos propagate in the loop as shown in the Fig|21 When the SUSY breaking mass 
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Figure 2: Chargino (a) and neutralino (b,c,d) loop diagrams to generate the LFV processes. 
The off-diagonal elements of rrij^ and m|, come from mj and m|, respectively. The mark • 
stands for a chirality flipping. The chirality flipping for Bino diagram (b) is given as m|^ = 
AeVd — fitanjS Me- There are diagrams in which the chirality is flipped in the external lines, 
but contributions from such diagrams are small. 

terms and the A-terms violate lepton flavor, the branching ratio of the LFV processes can 
be comparable to the experimental results. So the flavor structure of the SUSY breaking 
parameters is constrained PJ. 

When the SUSY breaking masses are universal and the A-term coefficient is proportional 
to the Yukawa coupling, there is no source for any LFV. However, even if there exists no 
LFV source in the SUSY breaking parameters at the cutoff scale, the sources for LFV can 
be generated through the neutrino Dirac Yukawa couplings as long as the coupling matrices 
are not proportional to the charged lepton Yukawa matrix. The Majorana couplings can also 
generate LFV. 

The RGEs above the scale of right-handed neutrino Majorana masses and and triplet Higgs 
fields are written in a proper notation as 

(4vr)^^m| = {Y^Y^ mj} + {Y^ , mj} + 3{hfl mj} (15) 

+ 2{Y,mlY^ + m\Y:Yl + A,A\) + 2{Y,m];Yl + m^jX + ^.4) 






+ 4(/r^(m^) Vr + fR'fR<- + ^/?^/^) 

-ft 

where {X, Y} = XY + YX and 5* is a trace of SUSY breaking masses with hypercharge weight. 
As it has been emphasized, the SUSY breaking scalar masses are universal due to the 
geometrical setup. On the other hand, if the F-terms of [/-moduli are zero, the A-terms 
are proportional to the Yukawa couplings. However, the non-zero values of F^ provide a 
source of LFV in the form as shown in the Eq.Q. Let us see the ^/-moduli contribution 
in the basis where the charged lepton Yukawa matrix is diagonal: UlYeU^ = Y^^^^. Since 



Ulx^ = (0, 0, Jaf + bl + cl)^, the derivative of the rank 1 part (Vo)ij = xfxf is written as 



/ X \ 



UliduYoW^^ = VI 



e\ 



re\ 



Vn' , (17) 



X 

y X X X y 

where the non-zero values in the entries are shown by x. Note that the (1,3) and (3,1) elements 
in the above matrix become zero when tan6'f^'^ = 1. Since the mixing angles in V^^ are 
small, the elements Aij {i,j < 2) can be small while the elements Ai^ and A^i can be large. 
In fact, the experimental bounds Br(/i — > 67) < 1.2 x 10^^^ [15j and the EDM of the electron 
\de\ < 1.6 X 10"^"^ 24] provide the most severe constraint on A12 and ImAn. Due to the 
structure of Eq.(fT7|). Ai^ and A^i from the U moduli contribution can be the sources of LFV 
while keeping the elements Au and An to be small. The [/-moduli contribution can generate 
the off-diagonal elements of SUSY breaking scalar mass squared matrices through the RGEs. 

In the minimal SUGRA, the non-proportionality of the A-term never develops. In the 
intersecting D-brane models, the non-minimality of the A-terms can be included when F^ 7^ 
while the SUSY breaking scalar masses have degeneracy for different generations at the 
cutoff scale due to the geometrical setup of D-branes. Due to the particular form of the U- 
moduli contribution as shown in the Eq. (fT7j) . the (1,3) element can be larger than the usual 
hierarchical assumptions for the non-minimal A-terms. The RGE effects are not decoupled till 
the electroweak scale and due to this the off-diagonal elements for both left- and right-handed 
slepton mass matrices are generated. The right-handed off-diagonal elements are always larger 
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than the left-handed elements due to a difference in the coefficients of the terms involving A^s 
in the RGEs. 

We enumerate the sources of LFV in the SUSY breaking scalar mass matrices at the weak 
scale in the mSUGRA model as follows: 

1. Neutrino Dirac Yukawa coupling [IB] 

The neutrino Dirac Yukawa coupling Yi^ can generate the off-diagonal elements of left- 
handed SUSY breaking slepton mass matrix mj . Hence, the chargino contribution can 
dominate in the LFV processes. The RGE effects are decoupled at the right-handed 
neutrino Majorana mass scale. 

2. Majorana coupling for left- and right-handed leptons 

The left-handed Majorana coupling fi is needed in type II seesaw. The right-handed 
Majorana coupling fn also participates in the light neutrino mass when B ~ L charge is 
gauged. The Majorana coupling //, and fn can generate the off-diagonal elements of both 
left- and right-handed slepton mass matrices, rn?^, m\. The RGE effects are decoupled at 
the ^L,R mass scale, and the right-handed neutrino Majorana mass scale for f^ coupling. 

3. SU(5) GUT 123 or the left-right unification ^ 

Since the right-handed selectron can be unified in the 10 dimensional representation of 
the SU{5) grand unification, the off-diagonal elements of right-handed selectron can be 
generated above the unified scale. The generated off-diagonal elements of m? are related 
to the CKM mixings. In the left-right unified models, the two Higgs bidoublets are 
needed to generate the CKM mixings, and the two different Yukawa matrices are sources 
of off-diagonal elements for both left- and right-handed sleptons. We do not discuss these 
sources in this paper. 

5 Numerical studies 

In this section, we will show the numerical calculations of the branching ratio of the LFV decays 
and the EDM of the electron in the context of intersecting D-brane models. 

We set up the parameters to show the numerical results as follows. The charged lepton mass 
matrix is given as rank 1 matrix plus small correction. The rank 1 matrix is given as Eq.(|H)) 
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in the basis where hght neutrino mass matrix is diagonaL In the minimal brane configuration, 
such as shown in the Fig.l, the parameters are given 



ai : bi : Ci = ^ 



e 




{t):§ 



-\ + e 




(t):^ 







(t). 



For the calculation, we use e = 0.1 and t = 1.5. Then 6^ = 47° and 9^ = 37°. For simplicity, 
the Yukawa matrix is assumed to be symmetric. Then [/-moduli contribution of the A-teim 
which is proportional to the derivative of Yukawa coupling is calculated in the basis where the 
charged lepton Yukawa matrix is diagonal 



A 



u 



cA, VI 



0.22 \ 
0.26 
V 0.22 0.26 0.84 j 






(19) 



where Aq is a dimensionful coupling coefficient and c is a coefficient. If F^ = 0, c = 0, the 
trilinear coupling is given as A^, = AqY^. + A]! . More precisely, the [/-moduli derivative of 
the correction matrix 6Y may also contribute, but we neglect its contribution here since its 
[/-moduli derivative does not appear to be large and has a model dependence. We will choose 
the mixing angles in VI = V^ as 6^2 — O-l^ ^23 ~ O-O^ and Of^ = 0.005. There can be 5 phases 
in the unitary matrix VI up to an overall phase in general, but for simplicity, we assume that 
there is no CP phase in VI . The neutrino mixing angles are given in the Eq. (jl2|) . In the choice 
above, Ue3 = 0.07 and 6'soi = 33°. 

The neutrino Dirac coupling can be written as Y^ = Vl'^Y^^^^V^ in the basis where the 
charged-lepton Yukawa coupling is diagonal. In a general scheme, the unitary matrix V^'^ is 
completely free. For example, Vl'^* is the MNSP matrix when type I seesaw is dominated and 
the right-handed neutrino Majorana mass matrix is proportional to identity matrix. However, 
in the present scheme of "almost rank 1 Yukawa matrices" , the unitary matrix Vl'^ is close to 
an identity matrix like the CKM matrix, and V^" ^ VI when the Dirac Yukawa coupling Y^, is 
hierarchical like up- type quark masses. We will use V^"^ = VI to express the numerical result. 
As we have already noted, we use 6^2 = 0.1, ^23 — 0-05 and 6' 



'13 



0.005. 



5.1 LFV decays 

We plot the branching ratios oi t -^ fij and /i — > 67 in Fig. El For the SUSY breaking 



parameters, we assume that mj 



m 



etj 



ttiq 1 and A^, 



AqY^ + A^l at the cutoff scale M* as 
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10-^" r 



T 



OQ 



10"' r 



10-" r 



— n— case 1 
— o— case 2 
—A— case 3 
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Figure 3: The branching ratios of the LFV decays are plotted. In the plot, slepton mass at 
cutoff scale is varied with 50 GeV steps. The detailed parameters we used are given in the text. 
Dashed lines are drawn for the current experimental bounds at 90%CL. The lines are plotted 
for the following cases: (Case 1) Y^,^ = Yt and c = 0, (Case 2) Y^^ = 0.1 Yt and c = 0.1, (Case 
3) Y^^ = Yt and c = 0.1. c is the coefficient given in the Eq. (fT!I|) . 

we have mentioned. In the intersecting D-brane models, the SUSY breaking scalar mass is not 
necessarily universal for different representation, though the flavor degeneracy is achieved. We 
assume that the left and right scalar masses to be same just for simplicity. The cutoff scale is 
related to the string scale and the volume of the extra dimensions. We choose that M* = 10^'^ 
GeV in the calculation. We take gaugino mass M1/2 = 500 GeV at M*, Aq = 500 GeV and 
Higgsino mass yU = 500 GeV (we choose the signature of fi to make the SUSY contribution of 
anomalous magnetic moment of muon ^27^ to be positive). The value of tan/3 which is the ratio 
of the vacuum expectation values for Higgs fields is taken to be tan/3 = 50. The amplitudes for 
the LFV decays are naively proportional to tan/?, thus the branching ratios are ex tan^ /?. In the 
Fig|21 we vary mo in 50 GeV steps and the maximal value is rriQ = 1250 GeV (which corresponds 
to the lightest stau of about 1 TeV). In order to show the results clearly, we assume that the 
LFV sources are only the Dirac neutrino Yukawa coupling and the [/-moduli contribution in 
Ae- We neglect the sources arising in the Majorana couplings fEq. ()14|) ) by assuming them to be 
small. In the plot, Dirac neutrino coupling is the only source in the case 1. We take the largest 
right-handed Majorana mass to be 10^^ GeV. In the case 2, the [/-moduli contribution is the 
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Figure 4: Plots for the ratio of branching ratios. The hnes are plotted for the following cases: 
(Case 1) n^ = Yt and c = 0, (Case 2) Y^^ = 0.1 Fj and c = 0.1, (Case 3) Y^^ = Yt and c = 0.1. 

dominant source of LFV. The case 3 has both sources. It is easy to see that the neutrino Dirac 
couplings makes the branching ratio Br(r -^ fi'j) large. This is because that these couplings 
generate the off-diagonal (2,3) element of the left-handed slepton mass matrix and contributes 
to the chargino diagram of r — > yU7. The flavor violation source arising from the neutrino Dirac 
couplings can contribute to the /i — ^ 67 decay since the (1,3) element is also generated, but this 
element is smaller than the (2,3) element. If we switch on the [/-moduli contribution, the (1,3) 
element can be comparable to the (2,3) element and thus, the [/-moduli contribution increases 
the /i — > 67 decay rate more than the r — > /i7 decay rate. The reason for the behavior of 
the lines being different (in the FigOI) when mo is smaller than 400 GeV is that the Bino-Bino 
diagram dominates rather than the chargino diagram due to the large left-right mixings of the 
slepton. The qualitative behaviors of the cases 1 and 2 are not very different even if we change 
the numerical parameters, but the case 3 depends much on the initial condition such as Aq 
and /i since there can be a slight cancellation among the diagrams. The right- and left-handed 
lepton decays do not have interference, and a huge cancellations among the diagram for the 
branching ratios happen hardly. 

The branching ratio for each decay mode depends on the initial conditions. However, as 
shown in the FigEl the ratio of the branching ratios can be a good prediction for different 
LFV sources. In the figure, we use the same initial conditions as before. The ratio of the 
branching ratio is almost determined by the mixing angles in Vl'^ for the case 1 and the ratio 
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of the (1,3) and the (2,3) elements in the A^ for the case 2. Therefore, if all the branching 
ratios for r -^ yU7, t ^ ej and fi —>■ ej are measured, we can obtain important information to 
identify the LFV source. In fact, the following relations are satisfied approximately: Br(/i — >■ 
e7)/Br(r -^ fi'y) ~ {Of.^)'^ /Bt{t -^ fiu^^Ur) for the case 1, ~ (Aei3Me33)^/Br(r -^ liv^Vr) for the 
case 2 and Br(r -^ e7)/Br(r -^ /i7) ~ (^13/^23)^ ^r the case 1, ~ {A^i-i / Ae2-iY ^^^ case 2. In 
the case 3, the LFV sources are mixed, we do not have such simple expressions. These ratios 
of the branching ratios do not depend much on the initial conditions such as mo, ^o^ ^^1/2, A'- 
and tan (3 if the chargino diagram provides the dominant contribution. When the sleptons are 
light and the left-right mixing becomes large, the Bino diagram can contribute to fi —>■ e'-f and 
bends the lines for the ratio Br(/i -^ e7)/Br(r —>■ fi'-f) in the FigEJfor smaller mo. 

The large Majorana couplings / can contribute to the LFV decays due to its off-diagonal 
terms. If the type II seesaw dominates the neutrino masses, the ratios of the branching ratios 
are almost determined by the neutrino mixings when the Majorana couplings are the dominant 
sources of the LFV decays. In this case, the ratio Br(/i —>■ e7)/Br(r —>■ fi'-f) is about U^^/Bt{t -^ 
HP^Ut) while the ratio Br(r -^ e7)/Br(r -^ fi'y) is about {Ues/Uf^sY ~ (^12)^- The first value is 
similar to the pure [/-moduli case (case 2) and the second value is similar to the Dirac neutrino 
case (case 1). Thus the observation of the ratios can sort out the LFV sources. 

Usually, the 13 mixing is smaller than the 23 mixing in V^^ or [/£, even if we use a different 
setup, and thus the r — >• 67 decay rate is expected to be smaller than the r — ;> /i7 decay rate. 
However, if the [/-moduli contribution dominates, those two decay rates can be comparable 
since the (1,3) and the (2,3) elements in A^! are comparable. So measuring the ratio Br(r -^ 
e7)/Br(r -^ fi'j) is very important to see the presence of the [/-moduli contribution. 

5.2 EDM 

The other important observables to select the sources of LFV are the EDMs of the electron and 
the muon. If the trilinear scalar coupling Aq and the Higgsino mass /i are complex parameters, 
the EDMs can be large even if we do not have any source of LFV violation. However, if those 
are complex in general, the EDM of the electron can be too large compared to the experimental 
bounds when the slepton masses are less than around 1 TeV J2H1- Thus a cancellation is needed 
in that case to satisfy the bound [221 • It is unnatural to have cancellations for both the electron 
and the muon and thus the muon EDM will be large enough to be detected in the future 
experiments |20J. It is often assumed that the Aq and /i are real to satisfy the experimental 
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Figure 5: Plots for electron EDM and branching ratio of /i — >^ 67. In the plot, the slepton 
mass at the cutoff scale is varied with 50 GeV steps. Dashed lines are drawn for the current 
experimental bounds at 90%CL. 

bound naturally. In this case, the amount of EDMs are related to the source of LFV. 

Let us suppose that Aq and /i are real and all the CP phases are in the Yukawa couplings. 
The EDMs are imaginary part of the amplitude of the loop diagram. Since in the diagram, 
where the chirality flipping vertex does not include CP phase in generation mixings, the Bino- 
Bino diagram dominates the EDM calculation. The imaginary parts of (Ae)ii and {Ae)i2ml2i 
etc in the basis where the charged-lepton matrix is real and diagonal can be also important. 
However, if the A-term is proportional to the Yukawa coupling, such imaginary parts are small. 
The electron EDM can be proportional to fitanf^rrir when the (1,3) mixings for both left- 
and right-handed slepton mass matrices are generated. If the neutrino Dirac Yukawa coupling 
is the only source of LFV, the off-diagonal elements of right-handed slepton mass matrix are 
very small and consequently, the EDM of the electron is small, de ~ 10^'^^ e cm. Even in the 
type I seesaw with generic right-handed Majorana mass matrix, the electron EDM is at most 
de ~ 10~^^ecm |31j. If the ^/-moduli is a source of LFV processes, the off-diagonal elements 
for both left- and right-handed slepton mass matrices can be generated and thus the electron 
EDM can be enhanced to reach the current experimental bound \de\ < 1.6 x 10~^^ecm |24j . 
Hence, the electron EDM is an important observable to see whether the LFV arises from only 
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the neutrino Dirac Yukawa coupling or not. 

In the intersecting D-brane models, if we assume that the F-term of moduli does not have 
any phase then Aq does not get any phase. The phase of the Higgsino mass depends on the 
model, but it may be related to the SUSY breaking parameters and thus can be real in such 
an assumption. On the other hand, the Yukawa couplings can be complex if we include the 
Wilson line phases in the theta function [9 . Then the [/-moduli contribution of the A term 
can be generically complex while Aq is real. 

We plot the electron EDM and the branching ratio of yU — > 67 in the Figl^l in the case 
where the [/-moduli contribution is large enough using the same input for the SUSY breaking 
parameters as before. In this case, the values different from what has be shown for the neutrino 
Dirac Yukawa couplings do not change the plot very much since the chargino diagram does 
not contribute to the EDM. In general each component of A^ can be complex independently. 
In the plot, we take the overall factor for the U moduli contribution c to be pure imaginary 
for simplicity. The EDM can easily saturate the current experimental bound. For larger mo, 
the AeVd contribution is larger compared to the fi tan (3 rrij. part in the left-right slepton mixing 
while for smaller rriQ (< 400 GeV), the /x tan /5 771,- contribution becomes larger than the A^Vd 
part. This is because /itan/JrriT- contribution needs triple mass insertion in the Bino diagram 
and thus the amplitude is suppressed by larger power of ttio than the A^Vd contribution which 
can be produced by a single mass insertion. 

The presence of large complex Majorana couplings can also saturate the EDM bound when 
H tan j3 iTLr contribution is large for light sleptons. The A^Vd part is small when [/-moduli contri- 
bution is absent. When sleptons are heavy (depending on /itan/?), the /itan/Jm,- contribution 
is suppressed due to the triple mass insertion and thus the EDM becomes smaller for a fixed 
Br(/i -^ 67), comparing to the case when the [/-moduli contribution is dominant. 

The muon EDM is rf^ ~ 10~^^ — 10~^^ecm as long as the bound for de is satisfied and 
there is no huge cancellation in de- The ratio of the EDMs does not depend on the ratio of the 
corresponding charged lepton masses, d^j_/de 7^ m^jme. 

6 Conclusions 

We have discussed the flavor sector in the intersecting D-brane models. In the D-brane models, 
the low energy effective action for the zero modes such as particles in MSSM can be calculated 
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using the geometrical parameters. The neutrino mixing is elegantly realized in the context of 
the almost rank 1 Yukawa matrix. The flavor degeneracy of SUSY breaking scalar masses are 
realized when the generation is simply replicated at the intersection of the D-branes. The non- 
proportional part of the scalar trilinear coupling is obtained when the F-term of the tZ-moduli 
is not zero. 

Emphasizing the flavor degeneracy of SUSY breaking scalar masses in these models, we 
study the lepton flavor violating processes. The [/-moduli contribution of the scalar trilinear 
coupling can be the source of lepton flavor violation as well as the neutrino Dirac and Majorana 
Yukawa couplings which are included in the MSSM plus right-handed neutrino. 

We calculate the branching ratios of the LFV decays and the EDM of the electron dg.- The 
observations of the d^ and the ratio of the branching ratios for /i -^ e'j, t —* fi'j and t ^ e'j 
are important to sort out the sources of LFV as shown in the FigHJ 

The bound for the EDM of the electron can be improved to de ~ 10~^^ ecm in the planned 
experiments |^. Under the assumption that the SUSY breaking parameters and the Higgsino 
mass are real, it is possible to see whether we have any source of LFV in addition to the 
Dirac neutrino Yukawa couplings in mSUGRA. The [/-moduli contribution in our scheme can 
saturate the current experimental bounds for d^ < 1.6 x 10^^'^ ecm as well as the Br(/i -^ 67). 
The current bound for the branching ratio is Br(|U —>■ 67) < 1.2 x 10"^^ and it can go down to 
~ 5 X 10~^^ in the near future [S3]- 

The [/-moduli contribution in the trilinear scalar coupling can make Br(r -^ e7)/Br(r -^ 
H'-f) to be order 1. In the models where the Dirac neutrino or the Majorana couplings are the 
primary sources of LFV, this ratio is much smaller than 1. In order to completely sort out the 
LFV sources from the ratio of the branching ratios, the r LFV decays with branching ratio 
10~^ — 10^^° need to be at least measured. At present, the upper bound on the branching ratio 
of r — > /i7 is 6.8 x 10~^ 15j and this bound can be improved to 10^^'^ at the ILC-Super B j34j . 

Acknowledgments 

We would like to thank Justin Albert for providing us the information on the future limit of 
the r — i> /i7 branching ratio at the ILC-Super B. 



18 



References 

[1] F. Gabbiani and A. Masiero, Nucl. Phys. B 322, 235 (1989); J. S. Hagelin, S. Kelley and 
T. Tanaka, Nucl. Phys. B 415, 293 (1994); F. Gabbiani, E. Gabrielli, A. Masiero and 
L. Silvestrini, Nucl. Phys. B 477, 321 (1996) |hep-p h/9604387|. 

[2] A. H. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982); R. Barbieri, 
S. Ferrara and C. A. Savoy, Phys. Lett. B 119, 343 (1982); L. J. Hall, J. D. Lykken and 
S. Weinberg, Phys. Rev. D 27, 2359 (1983); E. Cremmer, P. Fayet and L. Girardello, Phys. 
Lett. B 122, 41 (1983); N. Ohta, Prog. Theor. Phys. 70, 542 (1983). 

[3] M. Berkooz, M. R. Douglas and R. G. Leigh, Nucl. Phys. B 480, 265 (1996) 
|hep-th/9606T39i . H. Arfaei and M. M. Sheikh Jabbari, Phys. Lett. B 394, 288 (1997) 
|hep-th/ 9608167|. 

[4] R. Blumenhagen, L. Goerhch, B. Kors and D. Lust, JHEP 0010, 006 (2000) 
|hep-th/0007024| ; R. Blumenhagen, B. Kors and D. Lust, JHEP 0102, 030 (2001) 



|hep-th /0012156| ; G. Aldazabal, S. Franco, L. E. Ibanez, R. Rabadan and A. M. Uranga, 
JHEP 0102, 047 (2001) |hep-ph/001lT32| . 



[5] For a reivew, R. Blumenhagen, M. Cvetic, P. Langacker and G. Shiu, hep-th/0502005| 



^ 



R. Blumenhagen, B. Kors, D. Lust and T. Ott, Nucl. Phys. B 616, 3 (2001) 
|hep-th/0107'T38| ; M. Cvetic, G. Shiu and A. M. Uranga, Phys. Rev. Lett. 87, 201801 



(2001) |hep-th/0107143| ; Nucl. Phys. B 615, 3 (2001) |hep-th/010 7166|; D. Cremades, 



L. E. Ibanez and F. Marchesano, JHEP 0207, 022 (2002) |hep-th/0203160. ; M. Cvetic, 
P. Langacker and G. Shiu, Nucl. Phys. B 642, 139 (2002) |hep-th/0206TT 5|; G. Ho- 
necker and T. Ott, Phys. Rev. D 70, 126010 (2004) |hep-th/ 0404055) ; C. Kokorelis, 



hep-th/04062"58| |hep-th/0410134; Nucl. Phys. B 732, 341 (2006) |hep-th/0412035 



C. M. Chen, T. Li and D. V. Nanopoulos, Nucl. Phys. B 732, 224 (2006) |hep-th/0509059 



G. K. Leontaris and J. Rizos, Phys. Lett. B 632, 710 (2006) |hep-ph/0510230 



[7] C. Kokorelis, JHEP 0208, 018 (2002) |hep-th/0203187| ; J. R. Ellis, P. Kanti and 

D. V. Nanopoulos, Nucl. Phys. B 647, 235 (2002) |hep-th/0206087| ; M. Axenides, 

E. Floratos and C. Kokorehs, JHEP 0310, 006 (2003) |hep-th/0307255| ; C. M. Chen, 
G. V. Kraniotis, V. E. Mayes, D. V. Nanopoulos and J. W. Walker, Phys. Lett. B 611, 



19 



156 (2005) |hep-th/0501182| ; Phys. Lett. B 625, 96 (2005) |hep-th/0 507232 ; C. M. Chen, 
V. E. Mayes and D. V. Nanopoulos, Phys. Lett. B 633, 618 (2006) |hep-th/0511135 . 



[8] M. Cvetic, T. Li and T. Liu, Nucl. Phys. B 698, 163 (2004) |hep-th/040306ll ; F. March- 
esano and G. Shiu, Phys. Rev. D 71, 011701 (2005) |hep-th /0408059 ; JHEP 0411, 
041 (2004) |hep-th/0409T32| ; M. Cvetic and T. Liu, Phys. Lett. B 610, 122 (2005) 
|hep-th/ 0409032|; M. Cvetic, T. Li and T. Liu, Phys. Rev. D 71, 106008 (2005) 
|hep-th/ 0501041|; C. M. Chen, T. Li and D. V. Nanopoulos, NucL Phys. B 740, 79 (2006) 
|hep-th/0601064| . 



[9] D. Cremades, L. E. Ibanez and F. Marchesano, JHEP 0307, 038 (2003) |hep-th/0302"T05 



JHEP 0405, 079 (2004) |hep-th/0404229|. 



[10] M. Cvetic and I. Papadimitriou, Phys. Rev. D 68, 046001 (2003) |hep-th/ 0303*083 



S. A. Abel and A. W. Owen, Nucl. Phys. B 663, 197 (2003) |hep-th /0303127 ; Nucl. 
Phys. B 682, 183 (2004) |hep-th/0310257| ; D. Lust, P. Mayr, R. Richter and S. Stieberger, 
Nucl. Phys. B 696, 205 (2004) |hep-th/0404'l34l ; S. A. Abel and B. W. Schofield, JHEP 
0506, 072 (2005) |hep-th/0412206|; T. Higaki, N. Kitazawa, T. Kobayashi and K. j. Taka- 
hashi, Phys. Rev. D 72, 086003 (2005) plep^h/0504019j; M. Bertolini, M. Billo, A. Lerda, 



J. F. Morales and R. Russo, .hep-th/0512067| 

[11] B. Kors and P. Nath, Nucl. Phys. B 681, 77 (2004) | hep-th /0309T67|. 

[12] P. G. Caniara, L. E. Ibanez and A. M. Uranga, Nucl. Phys. B 689, 195 (2004) 
|hep-th/0311241| ; D. Lust, S. Reffert and S. Stieberger, Nucl. Phys. B 706, 3 (2005) 
|hep-th/ 0406092]; Nucl. Phys. B 727, 264 (2005) |hep-t h/0410074|; L. E. Ibanez, Phys. 
Rev. D 71, 055005 (2005) |hep-ph70408064i ; A. Font and L. E. Ibanez, JHEP 0503, 040 
(2005) hep-th/0412150|. 



[13] N. Chamoun, S. Khalil and E. Lashin, Phys. Rev. D 69, 095011 (2004) |hep-ph/0309"l69 



S. A. Abel, O. Lebedev and J. Santiago, Nucl. Phys. B 696, 141 (2004) |hep-ph/0312T57 



N. Kitazawa, Nucl. Phys. B 699, 124 (2004) |hep-th/ 0401096|; N. Kitazawa, T. Kobayashi, 
N. Maru and N. Okada, Eur. Phys. J. C 40, 579 (2005) |hep-th/0406TT5| ; T. Noguchi, 
hep-ph/0410302| 



[14] B. Dutta and Y. Mimura, Phys. Lett. B 633, 761 (2006) |hep-ph/0512TfT 

20 



[15] M. L. Brooks et al. [MEGA Collaboration], Pliys. Rev. Lett. 83, 1521 (1999) 
|hep-ex/990 5013 ; K. Abe et al. [Belle Collaboration], Pliys. Rev. Lett. 92, 171802 (2004) 
|hep-ex/031002 9 ; B. Aubert et al. [BABAR Collaboration], Pliys. Rev. Lett. 95, 041802 



(2005) |hep-ex/0502032j ; Pliys. Rev. Lett. 96, 041801 (2006) |h^x/0508012 . 

[16] F. Borzumati and A. Masiero, Pliys. Rev. Lett. 57, 961 (1986); J. Hisano, T. Moroi, 
K. Tobe, M. Yamaguchi and T. Yanagida, Pliys. Lett. B 357, 579 (1995) |liep-pli/950i407| . 

[17] J. C. Pati and A. Salam, Pliys. Rev. D 10, 275 (1974). 

[18] M. Cvetic, P. Langacker, T. j. Li and T. Liu, Nucl. Pliys. B 709, 241 (2005) 
|hep-th/0407178 . 



[19] M. Cvetic, P. Langacker and J. Wang, Pliys. Rev. D 68, 046002 (2003) |liep-th/0303208 |. 

[20] A. Brignole, L. E. Ibanez and C. Munoz, 'hep'-pli/ 9707209' 

[21] C. Giunti and M. Tanimoto, Pliys. Rev. D 66, 053013 (2002) |liep-ph/0207096| ; 
P. H. Frampton, S. T. Petcov and W. Rodejoliann, Nucl. Phys. B 687, 31 (2004) 
|hep-ph/0401206| ; Z. z. Xing, Pliys. Lett. B 618, 141 (2005) |hep-ph/0503200| . A. Datta, 



L. Everett and P. Ramond, Pliys. Lett. B 620, 42 (2005) liep-pli /0503222]; J. Harada, 
hep-ph/0512294| 



[22] P. Minkowski, Pliys. Lett. B 67, 421 (1977); T. Yanagida, in proc. of KEK workshop, 
eds. O. Sawada and S. Suganioto (Tsukuba, 1979); M. Gell-Mann, P. Ramond and R. 
Slansky, in Supergravity, eds. P. van Nieuwenhuizen and D. Z. Freedman (North-Holland, 
Amsterdam, 1979); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[23] J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980); R. N. Mohapatra and 
G. Senjanovic, Phys. Rev. D 23, 165 (1981); G. Lazarides, Q. Shafi and C. Wetterich, 
Nucl. Phys. B 181, 287 (1981). 

[24] B. C. Regan, E. D. Commins, C. J. Schmidt and D. DeMille, Phys. Rev. Lett. 88, 071805 
(2002). 

[25] L. J. Hall, V. A. Kostelecky and S. Raby, Nucl. Phys. B 267, 415 (1986); R. Barbieri 
and L. J. Hall, Phys. Lett. B 338, 212 (1994) |hep-ph/9408406j ; R. Barbieri, L. J. Hall 

21 



and A. Strumia, Nucl. Phys. B 445, 219 (1995) |hep-ph/9501334| ; J. Hisano, T. Moroi, 
K. Tobe and M. Yamaguchi, Phys. Lett. B 391, 341 (1997) |hep-ph/9605296 ; J. Hisano 
and D. Nomura, Phys. Rev. D 59, 116005 (1999) |hep-ph/9810479| . 

[26 



[27; 



[28 



[29 



[30 

[31 

[32 
[33 
[34 



K. S. Babu, B. Dutta and R. N. Mohapatra, Phys. Lett. B 458, 93 (1999) |hep-ph/ 9904366|; 
N. G. Deshpande, B. Dutta and E. Keith, Phys. Rev. D 54, 730 (1996) |hep-ph/95123 98|. 

H. N. Brown et al. [Muon g-2 Collaboration], Phys. Rev. Lett. 86, 2227 (2001) 
|hep-ex/01020T7| . 

E. Franco and M. L. Mangano, Phys. Lett. B 135, 445 (1984); P. Nath, Phys. Rev. Lett. 
66, 2565 (1991); Y. Kizukuri and N. Oshimo, Phys. Rev. D 46, 3025 (1992). 

T. Ibrahim and P. Nath, Phys. Rev. D 58, 111301 (1998) hep-ph/9807501| ; M. Brhlik, 
G. J. Good and G. L. Kane, Phys. Rev. D 59, 115004 (1999) |hep -ph/9810457| ; E. Acco- 
mando, R. Arnowitt and B. Dutta, Phys. Rev. D 61, 115003 (2000) |hep-ph/9907"446| . 



Y. K. Semertzidis et a/., hep-ph/0012087 



J. R. Ellis, J. Hisano, M. Raidal and Y. Shimizu, Phys. Rev. D 66, 115013 (2002) 
|hep-ph/0206lT0|. 



S. K. Lamoreaux, nucl-ex/0109014 



M. Grassi [MEG Collaboration], Nucl. Phys. Proc. Suppl. 149, 369 (2005). 

Talk by M. Giorgi at the 2nd Workshop on Super B-Factory, 16-18 March 2006, Fras- 
cati, Italy |http : //www . Inf . infn. it7c onf erence/superb06 /talks/Giorgi .ppt| ; 
The SuperKEKB can improve the branching ratio to several xlO^^ 
|http : //belle . kek . jp/superb/ 1 . 



22 



